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Abstract 

We analyze the tail behavior of the maximum N of {W^(£) — t 2 : t > 0}, where W is standard 
Brownian motion on [0,oo) and give an asymptotic expansion for ¥{N > x}, as x — > oo. This 
extends a first order result on the tail behavior, which can be deduced from Husler and 
Piterbarg (1999). We also point out the relation between certain results in Groeneboom 
(2010) and Janson, Louchard and Marhn-Lof (2010). 

1 Introduction 

The distribution function of the maximum of Brownian motion minus a parabola was studied in 
the two recent papers Janson, Louchard and Martin-Lof (2010) and Groeneboom (2010), both 
for one-sided and two-sided Brownian motion. The characterization of the distribution function 
is somewhat different in the two papers, but both characterizations (unavoidably) involve Airy 
functions. In this note we address the tail behavior of the distribution, a topic that was not 
addressed in these papers. 

The tail behavior of the maximum plays an important role in certain recent studies on the 
asymptotic distribution of tests for monotone hazards, based on integral- type statistics measuring 
the distance between the empirical cumulative hazard function and its greatest convex minorant, 
for example in Groeneboom and Jongbloed (2010). 

Let N be defined by 

N = max{W(t) -t 2 }, t > 0, (1.1) 

where W is standard Brownian motion on [0, oo). It can be deduced from Theorem 2.1 in 
Husler and Piterbarg (1999) that the distribution function Fn of N satisfies: 

1 f 8x 3 / 2 1 
l-F^-^expj- — j,*^oc. 

In section 2 we will give an asymptotic expansion of the left-hand side, which extends this result. 
The proof is based on an integral expression for the density, derived from Groeneboom (2010) 
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(which in turn relies on Groeneboom (1989)), and uses a saddle point method for the integral 
over a shifted path in the complex plane. As a side effect, it also leads to a clarification of the 
relation between the representations of the distribution, given in Janson, Louchard and Martin- 
Lof (2010) and Groeneboom (2010). 

2 Main results 

In the following, we will use Corollary 2.1 of Groeneboom (2010), which is stated below for ease 
of reference, specialized to the density of the maximum of W(t) - t 2 (instead of the more general 
W{t)-ct 2 ). 

Lemma 2.1 (Corollary 2.1 in Groeneboom (2010)) The density f of N is given by: 

f , ,/, x / ,- /« f°° Ai(e~ i7T / 6 u) Ai'(iu + 2 2 / 3 x) \) 
f N (x) = 2 2 > 3 <^ Ai(2 2 / 3 x) - 2 Re f e" m/6 y — i 1 du j ,00. (2.2) 

where Ai is the Airy function Ai, as defined in, e.g., Olver (2010) 1 . 

We deduce from this the following representation which is better suited for our purposes. 
Lemma 2.2 The density /jv of N is given by: 

2 2/3 / roc Ai(iu + 2 2 / 3 x) \ 1 f°° Ai(iu + 2 2 / 3 x) , 

f^ =— Re [l \ i(i ur d ")=¥^L \*«r ""■ x>0 - (2 ' 3) 

Proof. Integration by parts of the second term of (2.2) yields: 

f N (x) = 2 • 2 2 / 3 Re / AUiu + 2 2 / 3 x) — <^ — y — '- \ du , x > 0. 

\Jo v ' du y iAi{iu) j J 

Let the function h be defined by 

d j e~ i7r / 6 Ai {e-™/ Q u) 1 
~ du\ iAi(iu) J ' 

Using Olver (2010) 2 : 

Ai (V^/ 6 u) = Ai (e~ 2i7T / 3 iuj = ±e"^ /3 {Ai(tu) + iBi(iu)} , 

we obtain 

1 ( Ai' (iu) + iBi'(wt) (Ai(iu) + iBi(iu)) Ai'(m) 



h(u) 



2 \ iA\(iu) iA\(iu) 2 
and using the Wronskian Ai(z)Bi'(z) — Ai'(z)Bi(z) = 1/tt we conclude that 

1 



h(u) 



2-KAi(iu) 2 ' 

This gives the desired result. □ 



1 http://dlmf.nist.gav/9 
2 http://dlmf.nist.gov/9.2.Ell 
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Remark 2.1 Lemma 2.2 is in fact equivalent to relation (5.10) in Janson, Louchard and Martin- 
Lof (2010). The difference in the scaling constants is caused by the fact that they consider the 
maximum of W(t) — \t? instead of the maximum of W(t) — t 2 (see also section 3) and the fact 
that they integrate from — oo to oo (in that way also the imaginary part drops out). However, 
they arrive at this relation in a completely different way. So in this case we can go from Corollary 
2.1 in Groeneboom (2010) to the result in Janson, Louchard and Martin-Lof (2010), just by 
using integration by parts. This might serve as a first step in establishing the relation between the 
representations in the two papers. 

We are now ready to prove our main result. We will give two proofs, one based on the first 
equality in (2.3) and the other one based on the second equality. 

Theorem 2.1 Let N be defined by (1.1), and let fjy and Fn be the density and the distribution 
function of N, respectively. Then: 

ft) 

oo 

Mx) ~ 2 2 / 3 y^ J>A( X ), (2.4) 

k=0 

where 

$,(*) = — / e^- 22/3 *Hl- 3k dt. 

and dfc is given by: 

a = 1, oi = 2, a 2 = -20, a 3 = 560, a 4 = -25520, 
a 5 = 1601600, a 6 = -127568000, a 7 = 12287436800. 



(ii) More explicitly: 



fN{x) —exp\^-^-^ J 2^-^-, x^oo, (2.5) 



where the first coefficients are 

h — 1 h — !9,/q h — 3851 l _ 3380005 . /o l _ 6474441455 

b - 1, b 1 - 3 gV3, b 2 - -1536, o 3 - 221184 V3, b A - -- 

(Hi) 



14155776 



oo 

Cfc 



i tp r \ 1 J 8x3/2 I Ck 



where the first coefficients are: 



„ _ -I _ 19. /o „ _ _4535 _ 3869785 ./? . _ 7310315015 

M) — x 5 t-1 — 48 vo ' 02 1536 , t-3 221184 V °' 04 14155776 



Proof. Here we only derive the leading terms. Further terms in the asymptotic expansion are 
computed in the appendix. 
Let g(t) be defined by 



coo „—itu 
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By using the representation 

M{z) = — exp {it 3 - zt} dt. 

the second integral representation of /at(x) in (2.3) becomes 

f N (x) = 2" 1 /3 Re \ - / exp {it 3 - 2 2 / 3 xt) g(t) dt \ , (2.7) 

and this integral will be expanded for large values of x. For this we need an expansion of g(t) for 
large values of t. 

For large values of the argument \u\, and |arg(m)| < 7r, we have (see Olver (2010) 3 for the 
asymptotic behavior of the Airy functions): 

^i-^ ~ 2exp {l(iuf/ 2 } (1 + i)^Tu. (2.8) 

Hence, the dominant part of the integral in (2.6) is e^ u \ where 

<f>(u) = |(m) 3/2 - itu. (2.9) 
The derivative of this function vanishes at u s = —it 2 / 4, and we have at u s the expansion 

</>(u) = -±t 3 ~\{u- u s ) 2 + O ({u - u s ) 3 )/t 3 ) . (2.10) 

This suggests to take saddle point (for a changed integration road) for the integral in 

(2-6). 

We consider the following integration path: first the path Vi, going from to u s , and next the 
path V2, from u s to +00, into the valley of exp(|(iu) 3 / 2 ). That is, at +00 the phase of u is ^tt. 
See Figure 1, where we have shown the paths V\ and V2 for t = 1. 

We write g(t) = g\{t) + g2(t), where gj(t) is the contribution from the path Vj, j = 1, 2. Then: 

This shows that gi(t) is purely imaginary for t > 0. When we replace in (2.7) g(t) with gi(t), we 
see that this contribution to /n(x) becomes 

/3 



1 rcce"'>" I 

-2~ 1 / 3 Re I - / G(t) dt } , 

7T Jooe-Ki/a 



where G(t) is an analytic function which is real for positive values of t. The behavior of G(t) at 
infinity allows to take the contour along the imaginary axis. Integrating in this way, we obtain 



-2"V 3 Re jl jT° (G(iv) + G(-iv)) dv\ , 



and since G(iv) + G(—iv) is real for real v (because G(t) is a real function), this integral is purely 
real, and, hence, we can ignore this contribution. 



3 http://dlmf.nist.gov/9.7 
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Figure 1: Modification of the path of integration for the integral in (2.6) from (0, oo) to V\ U V2, 
for t = 1 with saddle point at u s = —i/4. 



Next we consider the integral over V2 ■ A parametrization of this path follows from the equation 
lm{(j)(u)} = lm{(f)(u s )}, and we see from (2.10) that lm{0(« s )} = 0. By writing u = rexp(i#), it 
follows that the path V2 can be described in terms of polar coordinates by 

3t 2 cos 2 6» „ , 

16sin 2 (^ + |7r)' 2^-^-6^ 

where for -> we have to apply l'Hopital's rule to obtain r = \t 2 . 

There is no need to follow this path for obtaining the asymptotic expansion of §2 (t) and for the 
sake of convenience we use the path from u s to u s + 00 parallel to the positive real axis. In this 
way we find the contribution from the path V2 , that is 

g— itu i>u s +co ^—itu 



^ = J V2 du = L du - (2 ' n) 

At the lower limit u s = —it 2 / 4 of the integration, we have by (2.8), 

exp{-|(m s ) 3 / 2 } 



irAi(iu s 

giving 



2(1 + i)V2u s = 2t, 



ru s +oo 

g 2 (t) ~ 2t / du 

J U s 



Next we neglect the O— term in the expansion in (2.10), and substitute u = u s + t 2 w, idGK. This 
yields 

Re(g(t)) = Re{g 2 {t)) ~ 2t 3 e~T2 t3 f°° e"' 3 ^ dw = t^V^y^, t -> 00. (2.12) 

Jo 



□ 
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Plugging the result in (2.12) into (2.7), and expanding at t = 2 4 / 3 y / x/3 yields: 

f 9-1/3 /-ooe™/ 3 ) 

Mx) „ Re | _ & eX p - 4'/3, t _ t 3} dt | 

( 9-1/3 /"ooe^/ 3 , . 

= Re J / t 3/2 exp f 1,3 _ 4 1/3 X A 



2 5/3 (x/3) 3/4 r ^3/2^1 ,00 

■exp^-— W exp|-2" 2/ VV3xj du 



Vtt I 3^/3 



f 8x 3 / 2 1 

= exp < — > , x — > 00. 

3 P \ 3^3]' 

This gives the leading term of the expansion; the further coefficients are computed in the appendix. 
Part (iii) follows by integrating this relation. 

(Second proof.) Again we only derive the leading term. We start with the second representation 
in (2.3). 

, 1 f°° Ai(iu + 2 2 / 3 x) 1 f°° Ai(u + 2 2 / 3 x) 

Mx) =¥7^L W - du =¥7HiL W x>0 - (2 - 13) 

and shift the integration path in the last integral to the path P along the line, parallel to the 
imaginary axis, and running from c — too to c + ioo, where c = |2 2//3 x. A similar path was used in 
the proof of (ii) of Corollary 3.4 in Groeneboom (1989). Asymptotically, as x — > 00, the exponent 
is now given by 

| (c + n) 3/2_ 2 (c + 2 2/3 x + n) 3/ 2) 

which equals -8x 3 / 2 /y/27 at u = 0. We get: 

1 [ c +™ M(u + 2 2 / 3 x) , 2c 1 /2^F e -8^/Vv^7 ,00 r 3u 2^3 ] 
MX) = 2VW C _ 100 Ai(«) 2 d "~ 2V3( C + 2^)1/4. L eXP \-^VVi))} du 

2 1 / 3 X 1 / 4 e -^ 3/ VV27 2 5/3 x l/4 4e -8o ; 3 / 2 /V27^ 



31/4 33/4 3 , ~->oo. 



□ 



Remark 2.2 In the first proof, the integral in (2.12) does not have a meaning for all complex 
values t. For example, we use values of t with phase tt/3 for which t 3 is negative. However, for 
all complex values of t with phases in [— tt/3, tt/3] we can give the integral in (2.12) a meaning 
by turning the path of integration in the w— plane such that the integral remains convergent. For 
example when arg(i) = tt/3 we can take arg(u>) = —tt/2. 

For two-sided Brownian motion we get similarly: 

Corollary 2.1 Let M be defined by 

M = max{VF(t) - t 2 }, 

where W is standard two-sided Brownian motion, originating from zero, and let fu an d Fm be the 
density and the distribution function of M, respectively. Then: 



6 



t f \ „ ML-M 8 ^ J 8x3/2 I 

f M {x) = 2f N (x)F N (x) — exp ( —£j=- } > x -»• oo. 

l-i^^-expj-^-j^oc. 

Proof. This follows from Corollary 2.2 of Groeneboom (2010), which gives the representation: 

f M (x) = 2f N (x)F N (x), x>0. 

□ 



3 Concluding remarks 

As pointed out to us by Svante Janson, the result implies certain facts for the moments of the 
distribution. For example, applying Theorem 4.5 in Janson and Chassaing (2004) together with 
Theorem 2.1 of the present paper gives: 

{EM r fl r ~ h$/2 r l 3 {r/e) 2 l 3 ,r^ oo. 
Theorem 2.1 can also easily be extended to a result for 

M c = max{W(t) - ct 2 }, 

by using the scaling relation: 

M c = c~ 1/3 M, 

see (1.7) of Janson, Louchard and Martin-Lof (2010). So, for example, Theorem 2.1 implies: 



W >, } = P {M>cV3 x }~-Lexp{-«} 



x — > oo. 



Also, by Lemma 2.2, the density of M c is given by: 

, x (4c) 1 / 3 / f°° AUiu + UcV^x) \ 
/M . (x) =U_ Be ^ Ai(au) 2 > *»),*><>■ (3-14) 

4 Appendix. Computing more coefficients of the asymptotic ex- 
pansions 

Recall, see (2.11), 

1 f e ~ itu 
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where we assume for the time being t > 0. We substitute u = t 2 v and write <j)(u) = —t 3 ip(v), where 
4>(u) is given in (2.9). This gives 



92(t) = - I 4 — dv, 

71 JPo c^O* v ) 



f 2 r e -t^(v 

b 2 e |(^) 3/2 Ai{it 2 v) 2 
where the modified path P 2 runs from —i/4 to — i/4 + 00, and 

We have iJ>(—i/4) = and substitute 



^(v) - -L = u, 2 . (4.15) 



14 



12 

Locally, we prescribe for this mapping 

w = ( v + i/4) + O ((v + i/4) 2 ) , v -i/4. 
Upon inverting we have -u = J2T=o^kW k , and the first coefficients are 

do = ~\h di = l, d 2 = \i, g? 3 = §, d4 = —§fi, d. 
The transformation (4.15) gives 

g 2 (t) = t 2 e-T2 t3 e- t3u,2 h(w) dw, (4.16) 
Jo 

where 

7ret^ 2 ^ 3/2 Ai(it 2 v) 2 

Because the argument of the Airy function is large for large values of t (for all values of w > 0) 
we expand the function h(w) first for large values of t. For this we need the well-known expansion 
(see Olver (2010) 4 ): 

Ai(^)~ - t- iTl E^^S^ Iph^l < vr, (4.17) 

where 

c = F /2 , 

and 

(2k + l)(2fc + 3)(2fe + 5) ■ ■ ■ (6fc - 1) T(3A; + \) 

Uk ~ 



(216) fc fc! 54 k klT(k + \Y 

The first coefficients are 

„ _ 1 _ _5_ _ 385 _ 85085 37182145 

' 1 72' " 2 — 10368' U3 ~~ 2239488' U4 ~~ 644972544' 

For /i(w) of (4.16) we obtain 

fc=0 s 



4 http://dlmf.nist.gov/9.7.E5 
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where the first coefficients are given by 

,„ _ 1 „„ _ 5 ... _ 155 „,, _ 17315 „„ _ 3924815 
Wq _ 1, Wl _ _ W2 _ W 3 - 0700075, w A - 



36' wz 2592' 279936' 40310784 ' 

We expand h(w) = YlT=o hk{t)w k and observe that the terms with odd index are purely imaginary 
and those with even index are real. We obtain 

Re(g 2 (t))~t 2 e-^ t3 J2h2k(t) / e^V* dw = \^ Vie'^ 3 £ { -%±h 2k (t). (4.18) 

Each Ii2k(t)/t is an infinite series in negative powers of t, representing an asymptotic expansion 
of these coefficients. For example, 

/*>(*) ~*(2 + fr 3 - ift-« + mpt- 9 + ...,) 

h 2 (t) ~ * (| - f^t" 3 + - ^|5450t- 9 + ...,) 

h (+\ + ( 140 1 36020 ^-3 6517750 ^-6 , 3140421550 +-9 , \ 
h 4 (t) ~ t (--^ + -gj— * 243— * + 2187 * + ■■■,)■ 

Using these expansions and rearranging the series, we obtain 



OO 



Re( 52 (t))~v^t 3 /V^ 3 ^^, (4.19) 

fc=0 



where 



a = 1, 01 = 2, a 2 = -20, a 3 = 560, a 4 = -25520, 
a 5 = 1601600, a 6 = -127568000, a 7 = 12287436800. 

We have derived the expansion in (4.19) with the assumption t > 0. As observed in Remark 2.2 
the integrals in (4.18) have a meaning for all phases of t € [— vr/3, tt/3]. Therefore, we can use the 
result (4.19) in (2.7) and obtain 

OO 

/ w (x)~2 2 / 3 ^^a,<D fc (x), (4.20) 

k=0 

where 

<i> fc(x) = _J_ / e it 3 -2V3^|_ 3fc ^ 

We continue the analysis considering the Airy-type integrals 

1 / I,„3_ 



A a (z) = — e* w - zw w a dw, (4.21) 

2t^ Joo e -Ti/3 

for large values of z. The contour should cross the real axis for positive values of w. We have 
* k (x) = (4/3) 5 / 6 -%/2-3fc (2 4 / 3 3- 1 / 3 x) , fc = 0, 1, 2, . . . . 
For nonnegative integer values of a the functions A a (z) are derivatives of the Airy function: 

A n (z) = (-l) n —Ai(z), n = 0,1,2,.... 
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For negative integer values of a they can be viewed as antiderivatives. We have, for example, 

/>oo 

A-i(z) = / Ai(t)dt. 

J z 

By using asymptotic expansions of the functions A a (z) (as given in the lemma below) and the 
corresponding expansions for the <3?fc in (4.20) we finally obtain after a few manipulations of series 

«*>~^«p(-f^)gjk- < 4 - 22 > 

where the first coefficients are 

f _ i f _ 19. /q f _ 3851 f _ 3380005 . A? f _ 6474441455 
/O - 1, 71— 4gVo, /2— T536, 73 - 2 21184 V <3 ' 4 _ 14155776 ' 

Lemma 4.1 The function A a (z) defined in (4.21) has the following asymptotic expansion: 



a/2— 1/4 — f 00 (°0 

A»W - * 9 7=^ E^T^ IpM<vt, (4.23) 



fc=0 

w/iere 

~ 3 

and, in terms of the Gauss hypergeometric function, 



c = h 3/2 , 



( a ) _ r(3/c + i; 



1-2^1 (-a,-2fe;i-3fc;3) . (4.24) 



fe 54 fc A;!r(fe + | 
Proof. Substitute in (4.21) w = \fzs and write the integral in the form 

i(a+l)„-C r ooeW3 



^a(«) = / e * W S a <fe 

Jooe-™/ 3 



where 

0( a ) = ±s 3 -s+§ = l(s + 2)(s-lf 



The substitution i = (s - l)y/{s + 2)/3) gives 

|(a+l) -C /**°° 

a*^) = — — / * ^(t) dt, m = s a 

2vr? J_ ioc 

By expanding ^(t) = EfcLo^fc**. 



fc=0 



The coefficients ^fc can be represented as a Cauchy integral: 
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where C is a small circle around the origin. By integrating in the s— plane, it follows that 

ds 



((s + 2)/3)( fc + 1 V 2 (s- 1)*+!' 

where C is a circle around s = 1 with radius less than 1. By using Olde Daalhuis (2010) 5 : 

(_i)kr(!fc + i-a) , , 
* = W r(| t+ |- a ) (-a, J* + i - a; -2) . 

This form of the Gauss hyper geometric function has to be modified when ^(k + 1) = a, but using 
Olde Daalhuis (2010) 6 , and taking into account that k is an integer, we can write 

, Fl (—a,—k; \k + \- a; -2) = J^J^-J , Fl (_, I - | t; 3) . 

This gives 

1 T(3A: + i) . 
Using the duplication formula (2/c)! = 2 2k k\ T(k + \)/^/tt we have proved the lemma. 



Remark 4.1 The Gauss hyper geometric function in (4.24) has a simple finite representation be- 
cause k is an integer. Also, for a = 1, 2, 3, . . ., giving results for the derivatives of the Airy function, 
we have simple relations. For example, when a = 1 we have uj^ = Uk(l + 6fc)/(l — 6fc) = — Wfc, and 
the ffe are used in the expansion of Ai'(z), see Olver (2010) 7 . 

Remark 4.2 Similar expansion as in (4.23) can be found in Drazin & Reid (1981) (pp. 465-478) 
for slightly different functions. 
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